/3-YbAlB 4 : a critical nodal metal 
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We propose a model for the intrinsic quantum criticality of /3-YbAlB4, in which a vortex in 
momentum space gives rise to a new type of Fermi surface singularity The unquenched angular 
momentum of the \ J = 7/2, mj = ±5/2) Yb 4/-states generates a momentum-space line defect in 
the hybridization between 4/ and conduction electrons, leading to a quasi-two dimensional Fermi 
surface with a kj_ dispersion and a singular density of states proportional to E~ x ^ 2 . We discuss 
the implications of this line-node in momentum space for our current understanding of quantum 
criticality and its interplay with topology. 



Since their discovery, heavy fermion materials have 
provided a wealth of insights into correlated electron 
physics. These materials contain a matrix of localized 
magnetic moments formed from /-electrons immersed 
within a host metal; at low temperatures the spin- 
quenching entanglement of the /-moments with the con- 
duction electrons gives rise to a diversity of low en- 
ergy ground- states, including anisotropic superconduc- 
tors, Kondo insulators and quasiparticles with effective 
masses 100s of times that of bare electrons [HE]. An im- 
portant sub-class of heavy fermion materials exhibit the 
phenomenon of quantum criticality, whereby the metal 
can be tuned via pressure, doping or magnetic field 
through a zero temperature second-order quantum phase 
transition, around which they show non-Fermi Liquid be- 
havior and predisposition to superconductivity [3H5]. 

The discovery [6 , 7 of an intrinsically quantum critical 
heavy fermion superconductor /3-YbAlE>4 has attracted 
much interest. /3-YbAlE>4 exhibits quantum criticality 
without tuning: in the absence of a field it exhibits non- 
Fermi liquid behavior, with a T 3 / 2 temperature depen- 
dence of the resistivity and a T -1 / 2 divergence in the 
magnetic susceptibility, but a magnetic field induces an 
immediate cross-over into a Fermi liquid (FL) with a T 2 
resistivity, in which the susceptibility diverges as B~ x l 2 . 
Intriguingly, T/B scaling of the free energy and magneti- 
zation has been observed over 4 decades in magnetic field 
[7], pin-pointing the critical magnetic field to within 0.1 
mT of zero while demonstrating that the Fermi temper- 
ature of the field-induced FL is the Zeeman energy [7 . 

In this paper, we show that the properties of this mate- 
rial can be understood in terms of a nodal hybridization 
model. In essence /3-YbAlE>4 resembles a Kondo insula- 
tor, but one in which the the hybridization gap between 
the conduction and f-electrons vanishes along a line in 
momentum space, producing a critical semi- metal with a 
singular density of states. 

There are three known examples of such nodal mate- 
rials: CeNiSn, CeRhSb and CeCu4Sn, in which the hy- 



bridization appears to vanish linearly along a line in mo- 
mentum space, closing the gap to form a heavy fermion 
semi- metal [8] . In /3-YbAlE>4 the unusual local seven- 
fold symmetry of the ytterbium (Yb) site surrounded by 
boron (B) atoms protects a "high-spin" \J = 7/2, mj = 
±5/2) state, in which the /-electrons carry a large un- 
quenched orbital angular momentum. The |5/2) state 
carries at least two units of unquenched orbital momen- 
tum orientated along the c-axis, yet plane waves carry 
no orbital angular momentum in the direction of mo- 
tion, so the /-state is protected from hybridization with 
conduction electrons traveling along the c-axis. This 
causes the hybridization to develop a singular structure, 
V(k) ~ (k x ±ik y ) 2 vanishing as the square of the trans- 
verse momentum with a double-vorticity associated 
with the two unquenched units of orbital angular momen- 
tum. The electrons and holes at the band-edge then form 
an emergent two-dimensional electron gas with a disper- 
sion proportional to the square of the hybridization, 

E(k) ~ |V(k)| 2 ~ (k ± )\ (1) 

giving rise to a quartic dispersion with a divergent density 
of states N(E) oc E~ x l 2 . It is the field-induced doping 
of this two-dimensional heavy band that accounts for the 
unusual field-tuned behavior in /3-YbAlB 4 . 

In /3-YbAlB4 the Yb atoms form a honeycomb lat- 
tice, sandwiched between layers of B atoms, with the 
Yb atoms sitting between a pair of 7- member B rings, 
giving rise to a local environment with local seven-fold 
symmetry [6], as shown in Fig. [I] We shall assume 
that the Yb ions are in a nominal Yb 3+ , 4/ 13 config- 
uration, with total angular momentum J = 7/2. Pho- 
toemission spectroscopy indicates a microscopic valence 
of 2.75 [9] due to moment- conserving valence fluctuations 
Yb 3 + ^Yb 2 + + e~. 

J =7/2 crystal field operators with 7-fold and time- 
reversal symmetries conserve total J z , splitting the 
J=7/2 Yb multiplet into four Kramers doublets, each 
with definite \mj\. The Curie constant and the 
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FIG. 1. Showing the seven- fold symmetric environment of the 
Yb 3+ ions (large spheres) in /3-YbAlE>4 , sandwiched between 
two heptagonal rings of B atoms (small spheres). The blue 
surface is the orbital distribution in the rrij = ±5/2 state. 

Ising anisotropy of the magnetic susceptibility of j3- 
YbAlE>4 are consistent with a pure Yb ground-state dou- 
blet | J = 7/2, mj = ±5/2) [10], a configuration that ex- 
hibits maximal hybridization with the seven-fold boron 
rings. This Ising ground-state is also consistent with the 
large anisotropic g-factor observed in electron spin reso- 
nance measurements on j3- Yb AIB4 [TTJ [12] . 

We model the low energy physics of /3-YbAlE>4 as a 
layered Anderson lattice [TO] , 

H = J2 CknC^CW ± J2 H m(j), (2) 

n,/c,<j j 

where the first term describes a tight-binding boron con- 
duction electron band where n is the band index, and 

H m {j) = Vo(c] a X 0a (j) ± h.c.) ± E f X aa (j), (3) 

describes the hybridization with the Yb ion at site j and 
the energy level Ef of the /-electrons. Here, X Qa = 
|4/ 14 )(4/ 13 ,a| is a Hubbard operator linking the 4/ 13 , 
mj = a = ±5/2 state of the Yb 3+ ion to the completely 
filled shell Yb 2+ state |4/ 14 ). The operator 

c 5«= E 4(Rjp) y~(r P ), (4) 

pe(l,14),cr 

creates a conduction electron in a Wannier state delocal- 
ized across the seven-fold boron rings directly above and 
below the Yb ion at site j, with local / symmetry and 
J z = a = ±5/2. The Rj p = Rj ± r p are the locations 
of the fourteen boron sites around the Yb site j (see Fig. 
1). The hybridization matrix, 

y aa (r) = chYtAr) = ^= (^If Jjj (r),(5) 
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where the C£ a = (3a — a, §|^,a) are Clebsch-Gordan 
coefficients for the Yb 3+ , a = ±5/2 configurations. 



We employ a slave boson decomposition of the Hub- 
bard operators, X 0a (j) = bjfj a , where bj and fj a are 
a slave boson and an Abrikosov pseudo-fermion respec- 
tively; in a mean-field approximation, 

H m (j)=V *[c}j ja + h.c] + E f f] a f ja + A (r 2 - 1)(6) 

where Vq is the quasiparticle hybridization, renormal- 
ized by the mean- field amplitude of the slave boson field, 
r = taken to be constant at each site. Ao imposes 
the mean- field constraint (rif) ± r 2 = 1, while the renor- 
malized position of the /-level Ef = Ao ± Ef. 

Next, we transform to momentum space and evaluate 
the form-factor of the seven- fold symmetric Yb-B cluster. 
To obtain a simplified model, let us assume a single band 
of dispersion ek hybridizing with the Yb atom. Rewriting 
the creation operator at a given boron site in terms of a 
plane wave state 4( R j'p) = (4A/")" 1/2 ^ k cj ca e" <k " R ^, 
and fj a = A/" -1 / 2 Xlk/ k " e2k Rj where J\f is the number 
of Yb sites, Eq. Q becomes: 

ct Q = (4A0- 1/2 E c L7.a(k)e- ik -^, 

kcr 

where the form-factor of the Yb-B cluster 

[7(k)U= J2 y^(r P )e~^. (7) 

p=l,14 

The mean-field Hamiltonian (6) can then be written in 
terms of the plane- wave and f\^ a operators as 

where all details of the hybridization are hidden in the 
matrix |V(k)] = |V r *7(k). Now in polar co-ordinates 

where we denote (cos (9, sin 6) = (cq,sq). The impor- 
tant point, is that the hybridization vanishes as sin 2 
along the c-axis. Now the effect of Fourier transform- 
ing in Eq. ([7]), is to replace the real-space argument by 
the momentum y(r) — >• 3^(k). To obtain an analytic ex- 
pression, we approximate the discrete sum over the po- 
sitions in the seven- fold B ring by a continuous integral: 
Y, p ?E± / If- We find that ^( k ) is proportional to 
a unitary matrix, 

v ^ = ^(-k<)> (10) 

where Vq — ^J~^ and 

ak = 6s'm(k z a/2)(k x ± ik y ) 2 J2(k±R) 

A = cos(k z a/2)(k x ± ik y ) 3 J 3 (k ± R), (11) 
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FIG. 2. (a) Showing dispersion around the c-axis, with an 
electron pocket at the T point and a hole pocket at the Z 
point, (b) Magnetic field fills the k\ band. 



where J n are Bessel functions of order n, R is the radius of 
the seven-fold rings and a is the distance between boron 
layers. Since J n (x) ex x n at small x, near the c-axis, the 
hybridization vanishes as k\, with a diagonal form 



V(k) 



(k x + ikyY 








As one proceeds around the c-axis, the phase of the 
hybridization advances by 47r, forming a double vortex 
in the hybridization along the c-axis. This vorticity is 
a consequence of angular momentum conservation about 
the c-axis: plane waves |k<j) traveling along the c-axis 
carry a spin angular momentum of ±| along the c-axis, 
and because the /-states are in an mj = =b|, angular mo- 
mentum conservation prevents the mixing of conduction 
and /-electron waves travelling along the c-axis. 

We can diagonalize the mean-field Hamiltonian, to ob- 
tain a hybridized dispersion 
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•£/)±[ 5 (c k - 



^ / ) 2 + |V(^)| 2 ]% (12) 



where \V(k)\ 2 = V§[\a^\ 2 + |/3 k | 2 ]. Fig. § illustrates the 
hybridized band-structure. Near the c-axis, the squared 
hybridization vanishes as V(k) 2 = A(k z )k\. The disper- 
sion in the vicinity of the c-axis is then given by 



region where e{k z ) < 0. In the case where e{k z ) changes 
sign along the c-axis, a two dimensional electron and hole 
band is formed above and below the /-level. 

To explain the intrinsic criticality of /3-YbAlE>4 we con- 
jecture that the /-level is pinned to zero energy Ef = 0. 
A heuristic argument for this assumption, is to regard 
/3-YbAlE>4 as a Kondo insulator in which the nodal hy- 
bridization closes the gap along the c-axis, pinching the 
/-level in the gap at precisely zero energy. At the cur- 
rent stage of understanding, this assumption is purely 
phenomenological, a point we return to later. 

If Ef = 0, the density of states for this dispersing sys- 
tem is then given by N*(E) = J2± N£(E)0(±E), where 



7V*(£) = 2 / k± 



dkj_ dk z 



1 



E\T± 



(13) 



where 
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0[=fe(k z )] determines the 



characteristic scales T for the electron (+) and hole 
(-) branch of the dispersion. Powerlaw scaling will ex- 
tend out to characteristic Kondo temperature Tk of the 
system, so that the total weight x of f-electrons con- 
tained within the divergent peak is 2x = j T ^ K N*(E) ~ 
4y/T K /T , giving T = 4T K /x 2 . 

If the /-level is pinned to zero energy, then at low tem- 
peratures a Fermi line of zero energy excitations forms 
along the c-axis. In a field, the Zeeman-splitting of the 
/-level induces a singular polarization of nodal electron 
and hole bands, broadening the Fermi line into a distinct 
tubular Fermi surface. When a field is introduced, a spin- 
polarized Fermi surface grows around the line-zero in the 
hybridization, giving rise to a density of states of or- 
der N*[^/hbB] ~ £? -1 / 2 , leading to a Pauli susceptibility 
that diverges as \ ~ B~ x l 2 . We call this field-induced 
Fermi surface transition a "vortex transition". Vortex 
transitions are reminiscent of a Lifshitz transition, but 
whereas Lifshitz transitions are point defects in momen- 
tum space [E2IT4], the vortex transition is a line defect. 

We can model the singular thermodynamics of the sys- 
tem with the Free energy 



/oo 
dEN(E) ln[l + e -P(E-9» B B a )} 



E(k ± ,k z )=E f 



-e{k z [ 



■7](k z )k[ 



where rj(k z ) = i^fj^ and we have assumed that \e(k z ) \ is 
large compared to |V(/c)|. In other words, the system de- 
velops an emergent two-dimensional Fermi surface, with 
a kj_ dispersion. A hole band is formed in the region 
where e(k z ) > 0, while an electron band is formed in the 



where 



*(») = 



— f 



dx 



E M 1 



-x — ya~\ 



=±5/2 



and T 1/2 = (1/2) ]T ± T± 1/2 . Fig. ^ compares the ex- 
perimental scaling curve [7] with that predicted by our 
simple model. However, while a qualitatively good fit 



FIG. 3. Theoretical fit (red line) to the measured field- 
dependent magnetization of /^-YbAlEU from [7] (gray dots) 
using formula [l4| with gmj — 2.85 and To = 6.65ey. 

to the observations is obtained using a gyromagnetic ra- 
tio consistent with the single ion properties of Yb in j3- 
YbAlE>4 , the characteristic energy scale required to fit 
the experimental results is To ~ 6.5eV, far greater than 
the characteristic Kondo temperature (~ 200K) of this 
system [7]. Using our relationship To = 4Tx/x 2 , we can 
understand this scale by assuming that about x ~ 0.1 of 
the f-spectral weight is contained within the vortex metal 
contribution to the density of states. 

We now turn to discuss some of the assumptions be- 
hind our model. One issue is whether the plane-wave de- 
scription of the vortex metal survives inclusion of band- 
structure effects. In this situation, angular momentum 
is only conserved modulo nft, where n is the order of the 
symmetry group of the Yb environment, requiring n > 5 
to avoid any admixture of \mj\ = 3/2,1/2 states into 
the perfect ±5/2 doublet. In a model of /3-YbAlE>4 , us- 
ing tight-binding coupling within the B planes and per- 
fect hept agonal Yb rings, the nodal structure does in- 
deed survive, as shown in Fig. [4j However more work 
is required to understand whether the nodes persist in a 
more realistic model of /3-YbAlE>4 . Another key assump- 
tion is that that the pinching of the hybridization gap by 
the node perfectly pins the /-level to the Fermi surface. 
Ultimately, this must arise from Coulomb screening, an 
effect that also needs inclusion in future work. 

Independent support for our phenomenological the- 
ory is provided by the locally isostructural polymorph 
a— YbAlE>4, which in contrast to the /3 phase, has a FL 
ground state [15] . Both systems display comparable char- 
acteristic "Kondo" scales Tk ~ 200K [9]. Recent exper- 
iments indicate in a magnetic field, a-YbAlB4 develops 
a two dimensional Fermi liquid at fields B > 3T [16], 
suggesting that a-YbAlB4 is a phase in which the /-level 
has become detached from the Fermi energy. 

Finally, we note that vortex structure in the hybridiza- 
tion suggests a kind of topological line defect in momen- 
tum space. In Kondo insulators, the hybridization van- 
ishes at the high symmetry points forming point defects 
[13], corresponding to a homotopy ^(H) = Z2. Vortices 



FIG. 4. Fermi surface from a tight-binding calculation. Note 
the cylindrical feature along the z-axis linked to an almost 2D 
surface in the k x — k y plane. For clarity, the 1st BZ has been 
shifted by ir/h to move the node into the center of the zone. 

in the hybridization suggest a further one dimensional 
homotopy, Ui(H) = Z. This is an interesting direction 
for future work. 
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